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Abstract 

Numerical simulations of Diffusion-Limited and Reaction-Limited Cluster- 
Cluster Aggregation processes of identical particles are performed in a two- 
dimensional box. It is shown that, for concentrations larger than a charac- 
teristic gel concentration, the morphology of the resulting spanning cluster at 
the gel time t g exhibits a crossover length L c between percolation (I > L c ) 
and aggregation (I < L c ). L c vanishes when increasing c, and, at a critical 
concentration value c p (where L c — > 0) the entire spanning cluster scales as 
the percolating cluster obtained by standard percolation. Even if for c > c p 
the long-range correlations are similar to that of percolation, the vanishing 
links in the structure suggest that an homogeneous regime appears at small 
scales. 

PACS numbers: 61.43.H, 64.60.A 
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I. INTRODUCTION 



The aggregation processes of small particles have attracted a great theoretical interest in 
the last decade due to their wide range of applications. On can distinguish Particle- 

Cluster aggregation, whose prototype is the Witten-Sander model ||, and Cluster-Cluster 
aggregation ||[7J . While the former model applies to field- induced growth processes, such as 
electrodeposition KJ, dielectric breakdown |TD[, viscous fingering |TT| , |T2"f , etc., the latter 
explains true aggregation processes such as polymerisation fl^fl!}] , oil in water emulsions [T3 



soot particles in flames flT^] , flocculation of colloidal particles ||TT|— pT|] . Among Cluster-Cluster 
aggregation processes one distinguishes diffusion-limited cluster-cluster aggregation (DLCA) 
H|7] and chemically-limited (also called reaction-limited) aggregation (RLCA) p2|-p^1 which 
correspond respectively to fast (fully screened) and slow (partially screened) aggregation of 
colloids |7| . 



The CCA processes leads to a flocculation regime when the concentration is smaller than 
a characteristic gel concentration c g , and to a gelation regime otherwise. In flocculation 
regime it remains a single fractal cluster at the end of the aggregation process with a fractal 
dimension D equal to 1.45 in DLCA case (and 1.65 in RLCA case). In the gelation regime 
it appears an "infinite" cluster at a given time t g (hereafter refered as the gelling cluster). 
When working with a finite (but large) box, such a cluster is usually defined as touching 
the box from edge to edge as in the percolation theory ||25|| . Since both physical problems 
are quite similar, it is tempting to explain the infinite cluster formation in DLCA by the 
percolation theory which excludes all dynamical phenomena and this has been done by 
several authors [^-|3l|. The fractal dimension of the infinite cluster of DLCA has been 
found to be different than the one of percolation [^B] , ^7| and it has been argued that this 
discrepancy is due to dynamical effects. However, recent aggregation experiments |32| on 
attractive coagulated particles suggest that percolation transition occurs at the particle 
concentration 0.42. Therefore the question is still open: are CCA models and percolation 
theory compatible? 
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In this paper we show that the infinite cluster obtained at t g in CCA models exhibits a 
crossover length L c between percolation and aggregation, for concentration values c larger 
than the gel concentration c g . Furthermore, we show that L c vanishes when increasing c, and 
for c ~ 0.5 the percolation regime shows up at all length scales in the infinite cluster. These 
results suggest that CCA processes can be viewed as an irreversible percolation phenomenom, 



as the invasion percolation model without trapping ||33| , |34[| and the model for a diffusion front 
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II. THE MODEL 

The two-dimensional Diffusion-limited Cluster- Cluster Aggregation model consists in a 
Monte Carlo algorithm which builds clusters on a lattice within a square box of edge length 
L. Initially particles are distributed randomly (but uniformily in the space) on the lattice 
sites up to volume fraction (concentration) c. In order to insure that the diffusion coefficient 
of the clusters varies as the inverse of their radius R, a particle (or a cluster) is chosen 
randomly according to a probability: 

n a 

where a (=—l/D) is the kinetical exponent. Then the cluster performs a translational 
motion by one unit (taking account Periodic Boundary Conditions, PBC) in any of the 
four directions ±1,±1 chosen at random. If the cluster does not collide with another, the 
displacement is performed and the algorithm goes on by choosing again another cluster. If 
a collision occurs between two clusters they stick together forming a new large cluster. In 
our simulation we have considered that collision occurs when a particle of a cluster tries to 
occupy a particle of another cluster. In that case the cluster is not displaced but a bond is 
established between the two contacting particles. If there are more than one collision at a 
given motion, only one bond is chosen at random. This trick has used by Kolb in a reversible 



diffusion-limited cluster aggregation model PB[. This important variant implies that there 
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are no loops (as in off-lattice DLCA ||36|| ) and that there is no intrinsic percolation at the 
begining of the process: the intitial concentration can be varied up to unity. 

For concentrations larger than c g there exists a gel time t g where a cluster becomes in- 
finite. This gelling cluster is stored for a numerical analysis, and to compare, we leave the 
aggregation process continuing up to the time where it remains only one cluster (hereafter 
refered as the final cluster). In the case of Reaction Limited Cluster-cluster Aggregation 
|22Vf24|, in addition to the algorithm described above, it is introduced a sticking probabilty 
p(-C 1). after a collision a new bond is created only if a random number (uniformly dis- 
tributed between zero and one) is smaller than p. Here, we have performed two-dimensional 
simulations in boxes of different sizes up to L = 240. 



III. RESULTS AND DISCUSSION 

Typical results that compare qualitatevely the resulting morphologies of the gelling clus- 
ter (G) and the final cluster (F) in DLCA are are illustrated in figure 1. In the case of the 
gelling cluster the other remaining clusters have been discarded and therefore are not shown. 
Note that the morphology of the cluster shown in figure la is strongly reminiscent of that 
of a percolating cluster as obtained in standard percolation theory f25[ when occupying ran- 
domly the sites of a square lattice with the percolation probability p c = 0.59273. To make 
quantitative comparaisons between G and F, we have calculated the mass M dependent box 
size I using the mass-counting algorithm |37|]. In figure 2 we show the log- log plot of M 



versus I for the same concentration than those considered in figure 1. The resulting slope 
indicates that M scales as l Dp and l d , for G (open circles) and F (black squares), respectively. 
D p (=1.89) is close to the fractal dimension of the percolating cluster in two dimensions, 
and d{=2) is the spatial dimension. 

If the mass of the gelling cluster G or the mass of the final cluster F scales as l D we 
might be tempted to conclude that the quantity m = should not depend on /. To 
test this we show, in figure 3, the curves giving m as a function of I. Figures 3a and 3c 
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corresponds to G for D = D p and D=d, respectively and figures 3b and 3d correspond to F 
for D = D p and D=d, respectively. One observes on these figures that m is independent on 
/ only in cases 3a and 3d and for sufficiently large concentrations, confirming our first guess 
of different fractal dimensions, 1.89 and 2.00, for G and F respectively. Figures 3b and 3c 
are shown here as conter-examples to illustrate the high precisions on these estimates for 
the fractal dimensions. In general, i.e. for not too large concentrations, one can define a 
c-dependent crossover length L c such that m becomes independent on I only for I > L c . For 
I < L c , another linear regime is observed, which is better extended at very low concentration, 
with a slope of about -0.4 and -0.5 in cases 3a and 3d, respectively, corresponding to the 
fractal dimension D ~ 1.89 — 0.4 ~ 2.00 — 0.5 ~ 1.5, close to the one of DLCA clusters. 
Therefore L c defines a change of scaling regime between DLCA for short lengths and either 
percolation (G) or homogeneity (F) for large lengths. In fact, L c should be proportional 
to the characteristic length correlation £ (or average size) of DLCA fractal aggregates ]36| . 



Note that in all cases L c vanishes for c ~ 0.5, and this is the reason of the absence of 
crossover in the curves reported in figure 2. 

When estimating the fractal dimension of the aggregates D agg in the range of lengths 
corresponding to the aggregation regime (7 < L c ), we found, in both cases G and F, that it 
increases significantly with concentration (the data are reported in figure 4). Those results 



confirm previous conclusions in two dimensions |38| as well as in three dimensions ||39|| . As 
a consequence the c-dependence of L c cannot be annalysed as a simple scaling relation. 
Moreover, the discrete values of / that require mass-counting calculations, and finite size 
effects, impedes us to determine L c within a sufficiently small range of error. Anyway, all 
the results depicted in figs. 2 and 3 suggest that, at least for I > L c the mass of the gelling 
cluster G scales with I as a standard percolation cluster. However such results are not 
sufficient to insure that their morphologies are the same. It is well known that there exist 
additional quantities (different and independent on the fractal dimension) to characterize 
a percolating cluster. In principle, the morphology is entirely characterized by an infinite 
set of exponents |25|;|4(|, but here we shall focuse on two particular exponents: the fractal 



dimensions of the "backbone" D bb j|T| and the fractal dimension of the "links" [ji2j| . The 
backbone of a cluster is the ensemble which remains after removing dead-ends (or dangling- 
ends). For a percolating cluster in two dimension the fractal dimension of the backbone 
D bb is 1.61. The links (also called red bonds) are the sites of the cluster that are singly 
connected, that is, if we take out a link the connection between the entire cluster is broken. 
For a percolating cluster in two dimensions the fractal dimension of links Di is equal to 0.75. 

In order to identify the backbone and the links of a given cluster, we have been obliged 
(due to the PBC considered in the CCA simulations) to span the cluster out of the box, 
taking account of PBC. Then, we have applied to the new cluster configuration (with size 
larger than the box size L, in almost all situations) a procedure suggested by Hermann et 
al. |^3| to identify the backbone and the links. Finally, the resulting backbone (and links) is 



unspanned and returned inside the original box. In figures 1, the sites depicted with strong 
grey and black colors denotes the backbone and the links, repectively. 

We have calculated the fractal dimension of the backbone D bb for both G and F clusters 
and different concentration values. The results are depicted in figure 5 for two box sizes 
length L =90 (black symbols) and 120 (open symbols). Note that for small concentration 
values the D bb values, for both G and F clusters, are quite similar. This can be understood, 
since, as mentioned above, the fractal dimension in the aggregation regime is almost the 
same in both cases. This conclusion is more convincing for L = 120, but for larger c values, 
D bb of the gelling cluster G becomes smaller than D bb of the final cluster F. Furthermore, 
in the G case, D bb reaches the value 1.61 (the same value than for a percolating cluster) 
for c ~ 0.5. For c > 0.5, D bb saturates to 1.89 (the value of the fractal dimension of a 
percolating cluster). In the F case, the same kind of saturation phenomenon occurs but, 
for c > 0.5, D bb saturates to the spatial dimension d. The increasing value of D bb suggests 
that the backbone structure is mainly reflecting the fractal aggregate structure (Which also 
increases with concentration as shown in figure 4). It is only when the crossover length 
L c vanishes that the backbone structure becomes characteristic of the one of a percolating 
cluster, since D bb varies from 1.61 to 1.89 as shown in figure 5. 



In figure 6 it is shown that for c smaller than about 0.5, the fractal dimension Di of the 
links is approximately equal to 1, suggesting that the mass of the links scales with I in a 
trivial manner. However, for c ~ 0.5, Di reaches the value 0.75 as in a percolating cluster. 
The Di undependence on the box size L suggests that there exist a threshold concentration, 
close to 0.5 where the entire G cluster scales exactly as a percolating cluster. When c is 
increased above this threshold, D\ vanishes, similarly to the percolation theory when the 
occupation probability p is increased above p c ||40|| . The vanishing links in the gelling cluster 
for large c values suggest that, at very small scales, the system becomes homogeneous, even 
if, at large distances, the system scales as a percolating cluster as shown in figure 3. 

In order to appreciate the degree of generality of our results, we have also performed 
some calculations in the reaction-limited case (RLCA). In figure 7 we have reported the 
m(l) curves and it can be shown that they exhibit the same qualitative behavior than in 
the DLCA case (figure 3). These results suggest that the percolation scaling could exist in 



other kinds of CCA processes such as the ballistic-limited |[44|| , the convection-limited fr4"5 
and the fluctuating bond K| aggregation models. 



IV. CONCLUSION 

In this paper we have shown that the infinite cluster obtained at the gel time t g in CCA 
models exhibits a crossover length L c between aggregation and percolation. Moreover, L c 
vanishes at a critical concentration c p (~ 0.5) where the mass of the entire system (and 
its backbone and links) scales as for a percolating cluster obtained at p c . For c > c p the 
percolation regime persists at least at large scales, because the vanishing links suggest that 
an homogeneous regime apperas at small scales. The value obtained here for c p is close 
to the value reported in the above mentioned experimental work |32| (~ 0.42), but typical 
statistical fluctuations to determine a critical value impedes us to insure if there are some 
relation between these two critical probabilities. It might also be worth finding a relation 
between c p and the critical probability of the site percolation threshold (p c = 0.59273) on 
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a square lattice. On the other hand, the fact that close to c p the fractal dimension of the 
links Di becomes equal to that of the percolating clusters suggests that the critical behavior 
in both cases are the same, since the correlation lenght exponent v is given by |^7[ . 
Preliminary calculations on a cubic lattice suggest that the results reported here are quite 
general and extend in three dimensions. 

One of us (A. H.) would like to acknowledge support from CONICIT (Venezuela). The 
numerical calculations where done on the computers of the CNUSC (Centre Universitaire 
Sud de Calcul), Montpellier, France, with support from CNRS. 
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FIGURES 

FIG. 1. Typical configurations for a (a) gelling cluster G and (b) final cluster F, for c=0.5 and 
L=120. The sites shown in dark grey and black belong to the backbone and the links, respectively. 

FIG. 2. Log-Log plot of M versus I for c=0.5 and L=240, for the gelling cluster G (open circles) 
and the final cluster F (black squares). These curves result from averages over 40 simulations. 

FIG. 3. Log-Log plot of m(= M{l)/l D ) versus I for L = 240 and different concentration values 
(0.1, 0.2, 0.3,..., 0.8, from top to bottom) in the DLCA case, (a) and (c) correspond to G-cluster 
for D = Dp and D = d, respectively, (b) and (d) correspond to F-cluster for D = D p and D = d, 
respectively. These curves result from averages over 40 simulations. 

FIG. 4. Fractal dimension of the aggregates D agg (I < L c ) versus c, for L = 90 (black squares), 
L = 120 (black triangles) and L = 240 (open diamonds), resulting from averages over 80, 60 and 
40 simulations, respectively. 

FIG. 5. Fractal dimension of the backbone versus c, for L = 90 (black symbols) and 
L = 120 (open symbols). Square symbols and circle symbols denote for the gelling cluster 
and the final cluster, respectively. These data results from averages over 80 and 60 simulations, 
respectively. 

FIG. 6. Fractal dimension of the links D\ (I < L c ) versus c, for L = 60 (open circles), L = 90 
(black squares) and L = 120 (open diamonds), resulting from averages over 120, 80 and 60 simu- 
lations, respectively. 

FIG. 7. Log-Log plot of m(= M(l)/cl D ) versus I for L = 120 and different concentration values 
(0.2,0.4,. ..,1, from top to bottom) for the RLCA case (with a sticking probability equal to 0.005). 
(a) and (c) correspond to G-cluster for D = D p and D = d, respectively, (b) and (d) correspond 
to F-cluster for D = D p and D = d, respectively. These curves result from averages over 20 
simulations. 
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